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ANNOTATION 


An  analysis  is  made  of  the  boundary  layer  effect 
on  the  electrical  characteristics  of  a  Faraday  MHD 
generator.  An  electrical  substitution  circuit  and  a 
method  for  calculating  its  elements  are  proposed.  An 
analysis  is  made  of  the  local  circuit  for  the  substi¬ 
tution  of  a  channel  in  a  gas-dynamic  approximation. 

It  is  noted  that  the  questions  of  the  flow  of  current 
through  a  turbulent  boundary  layer  require  careful 
examination.  It  has  been  shown  that  the  existing 
methods  for  the  calculation  of  real  resistance  of  a 
boundary  layer,  which  are  based  on  a  break  in  inte¬ 
gration  on  a  certain  arbitrary  boundary  from  the 
electrode,  give  results  which  are  strongly  dependent 
on  the  selection  of  this  boundary.  Therefore  it  is 
advantageous  to  conduct  the  performance  calculation 
of  an  MKD  generator  for  a  number  of  assigned  or 
experimental  values  of  equivalent  resistances  of  the 
substitution  circuit. 

The  real  resistance  of  boundary  layers  can 
considerably  reduce  the  output  power  of  a  small- 
scale  MHD  generators  (for  instance,  by  2-3  times). 

The  boundary  layer  effect  can  also  become  signif¬ 
icant  for  generators  of  large  dimensions. 

The  influence  of  near-electrode  resistance  on 
the  transverse  conductivity  in  the  channel  of  a 
Faraday  MHD  generator  when  the  Hall  effect  exists 
is  analyzed.  The  results  obtained  can  be  used  for 
the  evaluation  of  the  boundary  layeV  effect  on  the 
electrical  characteristics  of  an  MHD  generator  as 
well  as  for  the  construction  of  an  algorithm  for 
the  numerical  computation  of  flow  in  the  channel. 
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INTRODUCTION 

The  nulysi:;  of  the  expe?'iments  conducted  up  to  the  present  time 
with  MHF  ■  -  -.erators  on  combustion  prolucts  shows  that  the  output 
character.' "'.'cs  cf  tkv  r«nerators,  especially  small-scale,  are 
substantia  j.y  lower  than  calculated  from  the  one-dimensional  theory 
[1,  2].  Thu.  in  a  -rums-er  of  experiments  the  electric  power  output  is 
2-3  times  less  than  the  computed  value  [3],  Investigations  with  the 
help  of  probes  showed  that  the  distribution  of  electric  field  is 
nonur.Iform  in  interelectrode  interval  of  the  MHD  channel,  it  is 
possible  to  isolate  the  flow  core  and  the  near-electrode  areas  of 
change  in  potential  which  occurs  at  distances  of  the  order  of  the 
thickness  of  the  thermal  boundary  layer  on  the  electrodes.  Usually 
the  cathode  drop  in  potential  A  the  channel  is  of  the  same  order 

as  anode  drop  *va  ,  whereas  in  the  gas  discharge  aVc»M4  .  In 
connection  with  this  it  is  possible  to  assume  that  the  potential  drop 
in  the  channels  of  generators • i3  specified  by  the  boundary  layers  on 
the  electrodes  and  hydrodynamic  processes  play  a  substantial  role  here 
This  question  was  discussed  in  the  literature  and  investigated 
experimentally  [3,  9,  11-16,  19], 

In  experimental  generators  a  decrease  in  the  open-circuit  voltage 
has  also  been  detected  [11,  13,  1*0.  Since  the  induced  electric 
field  is  proportional  to  the  gas  velocity,  then  leakages  of  current 
can  be  expected  on  the  boundary  layers  on  the  insulating  walls  [10, 

*»  9  1 
jO  J  • 
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Experimentally  the  facts  described  above  show  that  surface  effects 
can  exert  a  noticeable  influence  on  the  electrical  characteristics  of 
generators  and  on  events  in  the  flow  core. 

Basically  this  influence  amounts  to  the  following. 

The  cooling  of  gas  in  the  boundary  layer  decreases  conductivity 
in  the  case  of  an  equilibrium  molecular  plasma  and  increases  the  real 
resxstance  of  the  boundary  layer.  On  the  other  hand,  this  same  effect 
facilitates  the  reduction  in  leakage  currents  on  the  boundary  layer  on 
the  insulating  wall. 


For  relatively  short  channels  whose  length  comprises  several  gages 
the  two-dimensional  effects  become  significant,  especially  when  the 
Hall  effect  exists,  [4,  6].  In  this  case  near-electrode  effects 
decrease  the  Hall  current  and  create  a  longitudinal  electric  field 
which  Interferes  with  effective  conductivity  in  the  flow  core. 

There  is  another  aspect  of  the  influence  of  surface  processes  on 
characteristics  of  MHD  channels.  Under  conditions  of  the  strong 
interaction  of  flow  with  magnetic  field  there  is  the  possibility  of 
boundary-layer  separation  from  the  electrode  wall  of  the  channel.  In 
this  case  an  internal  hydro-  and  electrodynamics  of  the  channel  are 
strongly  changed  [2,  22]. 

The  influence  of  near-electrode  effects  on  the  hydrodynamics  of 
the  channel  is  not  examined  here,  but  primary  attention  is  given  to 
the  calculations  of  near-electrode  resistance  during  the  diffuse  flow 
of  current  through  the  boundary  Jayer  and  leakages  along  the  insulating 
walls.  Furthermore  an  evaluation  is  made  of  the  influence  of  near¬ 
electrode  resistances  on  the  change  in  the  transverse  conductivity  in 
the  channel  of  a  Faraday  MHD  generator  with  continuous  electrodes. 
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SI.  THE  ANALYZED  MODEL  OP  THE  HETEROGENEITIES 
OP  PLOW  IN  AN  MHD  CHANNEL 

For  the  typical  conditions  characteristic  for  supersonic  equilib¬ 
rium  of  an  MHD  generator  the  Reynolds  numbers  per  unit  of  length 

Re*  *»•  iQ*  m'1  .  Since  the  length  of  a  gas-dynamic  circuit  (nozzle  and 
channel)  comprises  several  gages  (<  10),  then  the  boundary  layers  are 
apparently  turbulent,  and  the  channel  flow  -  undeveloped  (the  boundary 
layer  thickness  $  less  than  the  characteristic  transverse  dimension  of 
the  channel). 

Let  us  examine  the  diffuse  flow  of  current  in  the  boundary  layer. 

We  will  consider  that  pressure  is  constant  in  the  cross  section  of  the 

channel  and  velocity  distribution  and  temperature  correspond  to  a 

turbulent  boundary  layer  (a  two-layered  arrangement  of  boundary  layer 

was  examined).  Such  an  approximation,  which  subsequently  is  called  gas 

dynamic,  is  fulfilled,  if  the  joule  heating  and  pondermotive  forces 

are  small  as  compared  with  the  heat  flux  and  friction  on  the  wall. 

In  the  nucleus  of  boundary  layer  at  an  electric  current  density 
2 

'v  1  A/cm  and  small  Hall  parameters  this  approximation  is  apparently 
justified.  In  immediate  proximity  to  the  electrode  at  relatively  low 
wall  temperatures  it  is  necessary  to  consider  the  processes  of 
diffusion,  emission,  the  energy  balance  for  electrons  and  heavy 
particles,  and  others  (see  below)  [20,  21,  23]. 

1.1.  The  velocity  distribution  in  the  nucleus  of  a  turbulent 
boundary  layer  ir.  the  case  of  exponential  approximation  takes  the 
form  [8] 
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where  e=yA  ;  y  -  the  coordinate,  normal  to  the  wall;  ^  -  the 
thicKness  of  the  dynamic  boundary  layer.  The  exponent  IX  is  the 
function  of  the  Reynolds  number  and  far  from  the  separation  point  of 
the  boundary  layer  takes  values  in  the  interval  1/7-1/10.  In  a 
viscous  (laminar)  sublayer  the  velocity  distribution  is  linear: 


"•Tfcf9* 

where  u.-u^icp?  -  dynamic  velocity;  fa  -  the  local  friction  factor; 

iVrj  -  kinematic  viscosity  at  wall  temperature  Tw  .  The  thickness 
of  viscous  sublayer  fa  and  velocity  (]A  at  IJ  =  fa  are  determined  in 
the  following  manner 


°A  ujcj7?’  • 


1.2.  For  the  boundary  layer  at  the  Pranatl  number  Pr  4  i  far 
from  the  separation  point  the  temperature  field  is  connected  with  the 
field  of  velocities  by  the  following  correlation  [8]: 


t-Tu  It)"’  u 

To.-  Tw  " ' a  >  u. 


(4) 


where  A  -  the  thickness  of  the  thermal  boundary  layer;  ; 

Taw*  T»  1T^£>)  -  the  adiabatic  wall  temperature;  'fs  J  +  %  Q  71*  Ml,  ({fm)  > 
-  the  recovery  factor;  f=Cp/C¥  ~  the  adiabatic  exponent. 

Hence  the  distribution  of  static  temperature  in  the  boundary  layer  is 
expressed  thus: 


ggESff^-”^  — -r  .:v„ 


1.3.  The  transverse  conductivity  of  gas  0t  is  determined  by  scalar 
conductivity  (index  "O'1)  and  the  Hall  parameter  for  electrons  ^ 

b-lfy-  <6) 

In  the  boundary  layer  of  a  weakly  ionized  equilibrium  gas  £?0  and  0 
when  pressure  gradient  does  not  exist  depend  only  on  temperature  and 
in  the  case  of  a  constant  transport  cross-section  of  electron 
scattering  can  be  presented  in  the  following  form: 


0,(y)  *5,~IT/TJ*  */>[&(<-  ■  <7) 

m-fi-'J/T-)*  (8> 

where  Js  -  the  potential  (energy)  of  ionization  of  atoms  of  the 
addition. 

1.4.  The  appearance  of  heterogeneities  in  the  transverse  flow 
pattern  (velocity  and  conductivity)  strongly  complicate  the  problem 
of  performance  calculation  of  a  generator,  which  in  general  becomes 
three-dimensional.  Since  the  flow  is  undeveloped  it  proves  to  be 
advisable  to  separate  the  phenomena  in  the  flow  core  and  the  boundary 
layer  [2,  11,  22 j,  considering,  naturally,  that  the  flow  core  and 
boundary  layer  are  hydraulically  and  electrically  connected  with  each 
other. 


Usually  the  influence  of  the  boundary  layer  on  the  hydrodynamics 
of  the  channel  is  considered  by  the  reduction  in  the  through  section 
of  the  channel  by  the  magnitude  of  the  displacement  thickness 


and  by  the 


emergence  of  disturbances  in  the  flow 


core  with  boundary-layer  separation.  With  this  the  development  of  the 


boundary  layer  is  determined  from  the  change  in  parameters  in  the  flow 


core  along  the  length  of  the  channel  calculated  without  considering 


the  boundary  layer  effect  [22]. 
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To  account  for  the-  boundary  layer  effect  on  the  electrical 
characteristics  of  a  generator  various  schemes  have  been  suggested 
[9,  11,  12,  14-16].  As  a  rule  at  first  the  real  resistance  of  the 
boundary  layer  is  calculated  which  is  then  included  in  the  load 
resistance.  This  method  is  apparently  valid  for  channels  with  sub¬ 
divided  electrodes,  when  the  linear  deformation  in  heterogeneities 
witnin  the  limits  of  every  electrode  is  small.  For  the  determination 
of  the  influence  of  the  heterogeneities  of  flow  on  the  electrical 
characteristics  of  a  Faraday  generator  with  continuous  electrodes 
(  V,  »  Cons{4  )  under  conditions  of  a  powerful  change  along  the  length 
of  parameters  in  the  flow  core  as  well  as  in  boundary  layer  it  is 
expedient  to  examine  the  local  electrical  circuit  of  substitution,  and 
to  calculate  the  effective  characteristics  of  the  channel  by  means  of 
series  integration  along  the  length.  The  introduction  of  the  local 
circuit  of  substitution  makes  it  possible  to  take  into  account  the 
influence  of  the  heterogeneities  of  flow  on  parameters  in  its  nucleus, 
therefore  on  flow  as  a  whole. 

The  joint  calculation  of  the  flow  core  and  boundary  layer  should 
be  conducted  in  series  and  is  possible  in  practice  only  when  using  a 
computer,  whereupon  the  correlation  between  the  nucleus  and  boundary 
layer  is  considered  on  every  step  of  the  calculation  network. 

The  calculation  of  flow  in  the  flow  core  is  reduced  to  the  solution 
of  a  Cauchy  problem  for  a  system  of  ordinary  differential  equations. 

The  number  of  initial  conditions  include,  for  example,  the  load 
coefficients  for  the  nucleus  K*o(Xi)  or  voltage  on  the  boundaries  of 
the  flow  core  or  electrodes,  and  also  all  the  boundary  layer  character¬ 
istics.  As  a  result  of  numerical  solution  the  values  are  found  for 
all  variables  in  the  following  node  of  the  calculated  grid,  besides 

and  the  values  connected  with  K+,  ,  and  also  the  characteristics 
necessary  for  the  boundary-layer  calculation  (  ,  form-parameters 

y a $11  -g"—  ,  H=-JWV  and  others).  Further  the  numerical  solution 

of  the  boundary  layer  equations  is  carried  out  for  the  electrode  and 
Insulating  walls.  This  determination  of  the  boundary  layer  character¬ 
istics  makes  it  possible  to  calculate  (for  instance,  by  the  iterative 


method)  the  load  coefficient  /C„  which  satisfies  the  condition  of 
constancy  of  electrode  voltage  V}  =  V3i  =  const  .  The  complexity  of 
computation  of  is  connected  with  the  fact  that  the  near-electrode 
drop  in  voltage  in  general  is  the  function  of  current  density  and 
therefore  the  load  coefficient. 

After  the  agreement  of  electrode  voltage  the  calculation  procedure 
is  repeated. 

Below  the  local  circuit  of  substitution  is  analyzed  and  the 
evaluations  of  the  magnitudes  of  its  individual  elements  are  given. 
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§  2.  LOCAL  ELECTRICAL  CIRCUIT  POR  SUBSTITUTION 
OF  AN  MHD  CHANNEL 

2.1.  Let  us  examine  the  stationary  flow  of  a  viscous  compressed 
gas  (Px  4  <}  with  equilibrium  electrical  conductivity  in  the 
rectangular  channel  of  a  Faraday  MHD  generator  at  small  magnetic 
Reynolds  numbers.  We  will  consider  that  flow  in  the  flow  core  is  one¬ 
dimensional,  end  effects  do  not  influence  the  picture  of  current 
distribution  in  the  section  in  question,  boundary  layers  are  flat  and 
their  thicknesses  on  the  anode  and  cathode  are  equal,  and  the 
induction  of  magnetic  field  is  constant  in  the  section  of  channel. 

Figure  1  shows  the  picture  of  the.  -current  distribution  in  question 
in  the  cross-sectional  flow  of  the  channel,  and  Fig.  2  -  the 
distribution  of  velocity  and  temperature  in  the  boundary  layers..  For 
the  engineering  calculation  of  electrical  characteristics  of  thd  MHD 
channel  with  a  flow  which  is  heterogeneous  in  cross  section  the  local 
electrical  circuit  of  substitution  was  used.  The  precise  calculation 
of  the  distribution  of  fields  and  currents  in  the  cross  section  of 
the  channel  is  complex  and  requires  the  solutions  of  a  two-dimensional 
boundary  value  problem  [5,  7 ]. 

Figure  3  depicts  the  proposed  local  circuit  of  substitution, 
qualitatively  corresponding  to  the  picture  of  current  distribution 
which  was  shown  in  Fig.  1  (subsequently  all  values  are  related  to  a 
unit  length  of  channel,  for  example,  current  dl/dx  ,  resistance 
,  where  /  ,  %  -  the  integral  values  of  current  and 
resistance) . 
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.  Cross  section  of  MHD  generator  channel. 

(1)  electrode;  (2)  insulation  wall;  (3)  flow  core 


Pig.  3.  The  electrical  circuit  of  substitution. 


In  Figs.  1-3  the  following  designations  are  introduced:  &3  -  the 
thickness  of  the  thermal  boundary  layer  on  the  electrode;  fi  »  U  - 
interelectrode  distance  and  width  of  through  section  of  the  channel; 
&0  -  zero  line  (see  below);  £■  =  £m  -  U«,  B(/l~c  fii)  -  emf  generated 

in  the  flow  core;  &£  =  I  uBdy  -  emf  generated  in  the  boundary  layer 

on  the  electrode.  We  will  assume  that  in  the  area  +&})- 

-  (A/2-&3)  the  equipotentials  are  parallel  to  the  plane  of  the 
electrode  over  the  entire  channel  width,  and  the  leakages* of  current 
from  the  flow  core  can  be  significant  only  near  the  angle  of  the  . 
channel  under  conditions  of  operation  of  the  generator  close  to  no”-  ' 
load  running. 

Such  an  assumption  makes  it  possible  to  introduce  the  local  load 
coefficient  for  the  flow  core,  determined  analogously  with  the  one¬ 
dimensional  theory: 

Ki  mK„=Ei/U-t),  (9) 

The  circuit  of  substitution  (Fig.  3)  includes  the  resistors  of  the 


9 


flow  core  Ri  ,  load  RH  ,  near-electrode  resistors  R-a  >  Rc  and  three 
banks  of  resistors  which  shunt  the  load:  R^  ,  RK  ,  Ru  .  The 
resistors  of  group  Rsj  determine  the  leakage  of  current  7*  over  the 
boundary  layers  on  the  insulation  walls.  Effective  resistance 
R*  ~V>/73  ,  where  Vj  -  voltage  on  the  section  from  (-^4  )  to 

(  Vfc-4j  )  ;  Ru  determines  the  fraction  of  leakage  current  which 
does  not  reach  the  electrode,  and  Rti  -  the  fraction  of  leakage  current 
from  the  electrode  into  the  area  a/p.  -(a/p  -a  0)  • 

For  some  working  conditions  on  the  internal  surface  of  channel  a 
layer  of  electro-conductive  material  can  be  formed  (deposition  from 
the  flow,  products  of  the  reaction  of  the  material  of  the  electrodes 
with  the  gas,  etc.).  The  influence  of  such  a  layer  on  the  electrical 
characteristics  of  a  generator  can  be  formally  taken  into  account  by 
means  of  the  introduction  into  the  circuit  of  substitution  of  the  bank 
of  resistors  R*} ,  each  of  which  is  analogous  to  the  corresponding 
resistor  of  group  R»j  .  Resistance  Rv  is  the  resistance  of  leakage 
current  Ju  on  the  insulating  walls  of  the  channel.  Let  us  note  that 
resistances  Ra,  Rc,  Riz  cannot  be  determined  by  the  boundary  layer 
and  by  contact  resistance,  but  by  the  resistance  of  micro-arcs  which 
appear  in  a  number  of  cases  ac  large  electric  current  densities.  These 
resistances  can  depend  on  current  and  magnetic  field  considerably  more 
strongly  than  in  the  case  of  the  diffuse  flow  of  a  current  [18],  This 
case  is  not  examined  here. 

2.2.  The  line  in  the  boundary  layer  on  insulating  wall  far.- from 
the  electrode,  where  the  electric  current  density  is  equal  to  zero, 

n- 

let  us  call  the  zero  line.  Condition  JysO  corresponds  to  equality 
.  Hence  the  relative  distance  of  the  zero  line  from  the 
wall  is  determined  thus: 

A  o  _  u  tfa  Mn. 

Ki  *  <10) 

Since  the  strength  of  current  leakage  on  boundary  layer  is  proportional 
to  £©  ,  then  it  can  be  significant  only  at  >  and  A0  exceeds 

the  thickness  of  the  laminar  sublayer  $  only  at  0.8. 


i 


At  A0<  d'A  the  value  of  A0  is  linear  on  Ki  : 

'$*  ''n 


Ag  ~  ) 


(id 


In  the  mode  of  free  running  (  /?w**'oc>)  on  any  straight  line 
y- Const  should  satisfy  the  condition: 


«/t. 


\jH(Z)dz=0  or  /Zj./U, 


(12) 


i.e.  always  Ki^  i- 


The  density  of  transverse  current  is  determined  by  the  Ohm  law: 


jy(z)*6,(z}[U(z)B-Ey(z)l. 


(13) 


For  the  evaluation  of  the  open-circuit  voltage  we  will  assume  that 


“/*  .  ~  9*  , 

J  5,E»dz*Ey!  S.dz. 

0  ® 


In  this  case  the  open-circuit  voltage  is 

A/i 


k*  i  Wr  <*■ 


(14.) 


If  conductivity  is  constant  in  the  cross  section  of  the  channel, 
then  open-circuit  voltage,  with  other  conditions  being  equal,  is 
determined  by  the  velocity  profile.  The  decrease  in  conductivity  in 
the  boundary  layer  can  lead  only  to  an  increase  in  the  open-circuit 
voltage  as  compared  with  the  case  6  =  const  . 


In  the  case  of  the  exponential  approximation  of  the  velocity 
profile  with  index  n  the  expression  for  the  open-circuit  voltage  can 
be  written  as  (for  simplicity  it  is  considered  that  )  : 


si 

i 
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V kksU~B/l  {(i-  4-  t)+  It  nTT  (i-Tfh-cr)}.  (15) 

Figure  4  shows  the  change  in  the  dimensionless  open-circuit 

voltage  depending  on  the  relative  boundary  layer  thick 

ness  for  fl  =  1/2-1/10  at  CL  =  h.  . 


JLafA 


Fig.  1».  The  dimensionless  open-circuit  voltage 
depending  on  the  relative  boundary-layer  thickness 

BS/afa-comi,  a=h,  U • 

It  i3  evident  that  at  n  <  1/7  and  2  S/a  <  0.5  the  reduction  of  V** 
due  to  the  leakages  of  current  over  the  boundary  layers  on  insulating 
walls  is  less  or  of  an  order  of  15%.  Therefore  during  the  engineering 
calculation  of  the  nominal  mode  of  operation  of  a  generator  (  Ki^O.S  ) 
the  leakages  of  current  in  the  turbulent  boundary  layers  can  be 
disregarded. 


§3.  the  evaluation  op  the  influence  op 

VARIOUS  FACTORS  ON  THE  REAL  RESISTANCE  OP 
THE  BOUNDARY  LAYER 

The  calculation  of  near-electrode  resistance  represents  a  complex 
problem  and  the  solutions  found  relate  only  to  special  cases  [17,  20, 

21,  23].  In  the  wall  areas  of  the  boundary  layer  on  the  electrode 
the  electron  concentration  and  conductivity  can  be  decreased  strongly, 
which  leads  to  the  appearance  of  diffusion  flows  to  the  electrode  and 
an  increase  in  Joule  heat  release.  In  this  case  it  is  necessary  to 
solve  jointly  the  equations  of  diffusion  and  energy  for  electrons 
[20,  23].  At  distances  from  the  electrode  of  the  order  of  Debye 
radius  the  quasineutrality  of  the  plasma  is  disrupted  and  the  potential 
distribution  is  determined  by  the  Poisson  equation.  Let  us  note  that 
for  typical  operating  conditions  of  equilibrium  MHD  generators  the 
Debye  radius  at  T5  ~  1500°K  is  much  greater  than  the  mean  free  path 
of  electrons. 

In  engineering  practice  for  the  determination  of  near-electrode 
resistance  they  use  various  simplified  method,  the  applicability  of 
which  under  those  or  other  conditions  is  not  frequently  founded. 

During  the  analysis  of  near-electrode  resistance  usually  two  components 
are  separated  out,  one  of  which  is  specified  by  the  boundary  layer, 
and  the  other  -  by  contact  resistance  of  the  "plasma-electrode”  [9>  1*0* 
The  computation  of  the  resistance  of  the  boundary  layer  is  connected 
with  the  integration  of  specific  resistance  for  the  boundary  layer. 

In  this  case  in  a  number  of  works  the  lower  limit  of  integration  is 
assumed  equal  to  zero,  the  profile  of  reduced  temperature  is  usually 
taken  as  corresponding  to  the  exponential  velocity  profile  with  an 


index  equal  to  1.7,  and  the  influence  of  viscous  sublayer  and  changes 
of  Hall's  parameter  in  boundary  layer  are  not  considered  [11,  14] . 

Such  an  approach  can  lead  to  significant  errors,  especially  in  the 
case  of  "cold"  electrodes.  In  works  [9,  26]  on  the  basis  of  the 
comparison  of  calculation  with  experiment  (the  conclusions  of  the  works 
[21]  were  used)  it  is  recommended  to  take  a  lower  limit  during  integra¬ 
tion  over  the  boundary  layer  thickness  equal  to  the  Debye  radius.  In 
this  case  the  resistance  of  the  boundary  layer  should  not  depend  on 
the  density  of  current. 

Let  us  examine  the  influence  of  various  factors  on  the  magnitude 
of  the  real  resistance  of  a  turbulent  boundary  layer. 

Let  us  present  the  near-electrode  resistance  in  the  form  of  the 
sum  of  two  in  series  connected  resistances: 

Rj  S^3i 

where  R3t  -  the  resistance  of  boundary  layer  on  an  electrode  with 
thickness  )j  -  contact  resistance  of  the  "plasma  electrode, 

to  which  the  jump  in  potential  fl  V4  corresponds.  In  this  work  the 
calculation  of  R3Z  was  not  conducted,  but  the  potential  jumps  on  the 
cathode  A  Vcg  and  anode  aVat.  were  considered  known.  The  values  of 
and  &Vaz  can  be  calculated  during  the  clarification  of  the 
detailed  mechanism  for  processes  near  the  surface  of  the  electrode  or 
determined  from  experimental  findings. 

If  we  disregard  the  boundary  layer  effect  on  the  insulator  on  near- 
electrode  resistance,  then  R3t  can  be  expressed  in  the  following  form:1 

‘Here  the  distance  along  axis  "  y  "  is  counted  off  from  the 
surface  of  the  electrode. 


*■ 


f„°Se-‘(U><lS>  - 

*•/&» 


(20) 


wnero  £a  -  y/A3  . 


Some  authors  assume  [9,  21]  that  the  value  of  ?*  is  determined  by 
the  characteristic  value  of  the  order  of  the  length  of  Debye  shielding 
in  plasma  Zp  ' 


ZM=  !r(  Zb )/4  fte*nc  (tm)  * 


(21) 


However,  already  at  y  >zn  even  in  molecular  gases  there  is  the 
possibility  of  a  significant  breakaway  of  the  temperature  of  electrons 
Tfl  from  the  temperature  of  heavy  particles  under  the  effect  of  the 
electric  field  increasing  toward  the  electrode  (  £*({/)=:  jy/0(o)  in  the 

model  in  question  jt,  ~  const) .  We  evaluate  the  influence  of  these 

approximations  or.  near-electrode  resistance  Jjt  and  the  conditions 

under  which  the  model  approximation  Z9~ZB  becomes  unacceptable  due 
to  the  appearance  of  a  noticeable  breakaway  of  the  temperature  of 
electrons . 


3.1.  Let  us  examine  the  case  ?.9~lD  (we  consider  tha»t  ) 

Usually  Xb  and  it  is  possible  to  isolate  the  real  resistance’ of 
the  turbulent  nucleus  and  viscous  sublayer: 

Jst  ~fr  +? a-  (22) 

The  results  of  calculations  Q  for  the  typical  combustion  products 

of  hydrocarbon  fuel  [3,  22 j  (  ^  =  1.12;  =  1.5;  7^  *  3000°K; 

j9 4,  1)  for  7?  =  300°K-2100°K  and  ft  =  1/5-1/10  are  presented  in  Pigs. 

5,  6.  At  ft  1/7  resistance  fs,  decreases  by  an  order  with  an 

increase  of  TV  from  300°K  to  1800°K  (Pig.  5).  The  small  change  in  the 
exponent  fl  at  T'-'C  1000°K  exerts  a  noticeable  influence  on  the 
resistance  of  the  boundary  layer. 
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boundary  layer  on  the  temoerature  of  an  electrode 
(U  ~  V"  ). 


Designation:  [om.m  =>  ohm.m]. 


Figure  6  shows  the  influence  of  the  viscous  sublayer  on  while 
maintaining  a  similitude  of  profiles  of  the  given  velocity  and 
temperature  for  H  =  1/7  in  the  nucleus.  It  is  evident  that  simple 
exponential  approximation  of  the  velocity  profile  up  to  the  wall, 
usually  utilized  for  computation  [9,  14],  at  low  temperatures  can  lead 
to  an  understating  of  the  magnitude  of  fw  by  several  times  in 
comparison  with  the  two-layered  arrangement  of  the  boundary  layer  even 
at  small  0.005). 

3.2.  The  distance  from  the  electrode  at  which  breakaway  7*  becomes 
significant  can  be  evaluated  by  the  equation  of  energy  balance  of  the 
electron  gas: 

(23) 


If  the  effective  values  of  the  collision  rate  V«  »  the  coefficient  of 
inelastic  losses  and  the  atomic  masses  (of  the  molecule)  are 

known,  then  the  electric  field  in  which  the  assigned  breakaway  of  the 
temperature  of  electrons  (  Te-T)  appears  is  calculated  as: 


Expression  (24)  makes  it  possible  in  the  case  of  known  dependence  6(<J) 
to  jy  -Jyoo  determine  the  characteristic  distance  t*  from  the  electrode 
at  which  at  a  distance  </aZt  the  "heating"  field 
guarantee  the  assigned  breakaway  of  the  temperature  of  electrons.. . 

During  the  calculation  of  the  breakaway  of  electron  temperature  it  is 
assumed  that  the  diffusion  length  of  the  ions  ii.  [24]  is  less  than  the 
value  7e> .  Figures  7,  ?,  and  9  depicts  the  comparative  calculations 
for  a  two-layered  arrangement  of  the  boundary  layer  at  Z9  s ta  and 
at  (T'-Tt&vw'  w'  and  various  current  densities. 

At  Z?-Zt  )  the  value  of  Zt  }  Tw  -  assigned)  is 

determined  by  the  conductivity  Q  (Zt)  »  which  is  attained  at  an 
identical  temperature  at  different  distances  from  wall,  which  gives 
rise  to  a  weak  dependence  of  P»*  on  the  structure  of  boundary  layer. 


300  600  900  (900  1600  (300  9(00 

Pig.  6.  The  influence  of  the  viscous  sublayer  on  the 
resistance  of  the  boundary  layer  =  0.003). 


IS 


Fig.  9.  The  projection  of  the  line  of 
intersection  of  surfaces  *•  W/ty and 

*e‘/'(l i.T,)  on  th<!  Plane  |  ,  T,  (curve  1  - 

for  the  exponential  velocity  profile, 
curve  2  -  for-  the  two-layered  circuit. 


22 


I 


fc 


The  influence  of  the  temperature  profile  on  the  magnitude  of  Pfi  ^ 
similar  to  the  case  Zt-Za  . 

Therefore  during  the  calculation  of  J3^  one  ought  to  select  as  the 
lower  limit  of  integration  in  expression  (20)  the  maximum  value  from 
ZB  and  Zt  .  Since  the  influence  of  the  resistance  of  the  boundary 
layer  on  the  characteristics  of  the  generator  (see  also  §  5)  is 
especially  significant  only  at  I500°K  ,  then  in  specific  interesting 

cases  (  p  g  10  ata;  1-5  a/cm2;  fappm  20-50;  Zf>ta 

(Pig.  9).  Let  us  note  that  at  Z3=Ze  the  value  of  f3t  drops  with  an 
increase  of  current  density,  other  conditions  being  equal.  However, 
this  does  not  mean  that  the  complete  near-electrode  resistance  pa 
decreases,  i.e.,  resistance  Psi  also  depends  on  the  current  density. 


I 

'c 

i 
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The  analysis  conducted  shows  that  the  existing  methods  for  the 
calculation  of  the  real  resistance  of  the  boundary  layer,  which  are 
based  on  a  break  in  integration  at  a  certain  arbitrary  distance  from 
the  electrode  t;  ,  give  results  which  are  strongly  dependent  on  the 
selection  of  this  boundary.  Integration  in  formula  (20)  up  to  maximum 
value  from  Za  and  Ze  apparently  gives  an  evaluation  of  the  minimum 
value  of  the  r.ear-electrode  resistance.  Therefore  the  performance 
calculation  of  an  MHD  generator  is  advantageously  conducted  for  a 
series  of  the  assigned  (parametric  analysis)  or  experimental  values 
of  the  equivalent  resistances  of  the  circuit  of  substitution. 
Specifically  parametric  analysis  will  make  it  possible  to  determine 
the  permissible  values  of  internal  resistance  of  the  MHD  channel. 
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!  §  4 .  THE  INFLUENCE  OF  NEAR-ELECTRODE  RESISTANCE 

ON  THE  EFFECTIVE  TRANSVERSE  CONDUCTIVITY  OF  A 
FARADAY  GENERATOR 

It  is  known  that  in  real  MHD  channels  with  continuous  electrodes 
there  are  always  both  significant  near-electrode  resistance  and  two- 
dimensional  effects  in  the  distribution  of  potential  and  current, 
especially  when  the  Hall  effect  exists.  These  effects  can  lead  to 
the  appearance  of  a  longitudinal  (Hall)  component  of  electric  field 
and,  consequently,  to  an  increase  in  effective  transverse  conductivity 
in  the  flow  core.  The  influence  of  the  two-dimensional  nature  on  the 
characteristics  of  generators  with  a  short  channel  (number  of  gages 
K-L/ix  ^3 )  was  examined  in  a  number  of  works  [25»  4,  6].  However,  in 
this  case  the  influence  of  the  heterogeneities  of  flow  connected  with 
the  presence  of  boundary  layers  was  not  analyzed.  Near-electrode 
resistances  prevent  the  complete  closing  of  Hall  current  in  a 
generator  with  continuous  electrodes  and  thereby  lead  to  the  appearance 
of  a  longitudinal  electric  field  and,  consequently,  to  an  increase  in 
transverse  conductivity  in  the  flow  core  as  compared  with  the  case  of 
uniform  distribution. 

We  will  evaluate  the  influence  of  this  effect  on  the  characteristics 
of  an  MHD  channel  with  continuous  electrodes. 

Let  us  write  the  components  of  current  densities  in  the  following 
form: 


2^ 


(25) 


Combining  expressions  (26)  and  (27) »  we  obtain  for  relation  -jfa 


the  following  formula: 


The  expression  for  the  density  of  transverse  current  we  present  in 
the  form  of 


h=U„8U-K„)&>]l+  7^3^ 


This  expression  can  be  reduced  to  the  common  form: 


where 


6ip<p  5,  (i*  773*17x11 )  ~ 


the  effective  transverse  conductivity  of  the  channel. 


The  appearance  of  a  longitudinal  electric  field  leads  to  an  increase 
in  effective  conductivity  as  compared  with  its  value  from  the  cine-; 

dimensional  theory  (5^  for  value  ( / 1  j r'» )  •  Since  near-  . 

electrode  resistances  lead,  on  the  one  hand,  to  obvious  added  losses 
in  total  generated  power,  and  on  the  other  (as  also  common  end  effects 
[25,  6])  -  to  an  increase  in  effective  conductivity  in  the  flow 

core,  then  under  specific  conditions  there  is  the  possibility  of 
existence  of  optimum  in  the  dependence  of  the  complete  generated 
power  on  the  characteristics  of  the  near-electrode  effects. 
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Output  power  is  maximum  under  the  condition  of  agreement  of 
resistances  of  the  channel  and  load: 


where 


Ri 


*p<p 


With  a  small  change  in  conductivity  and  velocity  in  the  section  of  the 
electrode  with  the  length  ^ *  the  expression  for  maximum  power  takes 
the  following  form: 


Let  us  rewrite  this  formula,  using  the  correlation 


(33) 


P*PoV(<X,<t',P,b*), 


(3*) 


where 

PQ  *  UiBi 'a A L*t 


(35) 


The  condition  of  the  extreme  of  function  ip  on  parameter  d(dif/dd-0) 
leads  to  the  following  equation: 


o(*[2(f+P*+°(*  %;KP*]  +  4<*otK(1+p*)* 


(1+p*+  ofn  \  K)+2  0. 
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(36) 


V.  .  gttijjl 


Analysis  shows  that  this  equation  has  a  positive  root  in  the  case  cf 
fulfillment  of  the  following  conditions: 

'  ?*-&*•  <»> 

The  values  of  parameters  diK  and  ^  can  be  evaluated  by  means  of  a 
comparison  with  the  results  of  two-dimensional  calculation.  The 
computation  of  power  using  formula  (3*0  leads  to  satisfactory  agreement 

with  the  results  of  work  [4]  at  l/fo  ~p/&  for  3  =  0-10.  In  this  case 

inequalities  (37)  are  fulfilled  at 


i.e.,  for  /C£  ?.  at  p  2  »  <**  ^  q. 

Figure  11  shows  the  dependence  of  relative  power  Lfs.  p//jo  from  3 
at  l^fi^/4  '  K=5 ~;d=Z;4  .  Curve  1  corresponds  to  calculations  in 
the  one-dimensional  theory  (  c(  *  *  °°» .  ;  curve  2  - 

using  formulas  (35),  (36)  at  ;  curve  3  -  at  {(,»  2  •  It  is 

evident  from  this  figure  that  at  large  3  the  difference  from  the  one¬ 
dimensional  theory  is  considerable  and  at  3  *  3-2  the  increase  in 
output  power  as  compared  with  one-dimensional  calculation  comprises 
-1.5. 

Figure  12  depicts  the  dependence  of  on  ct  at  fa  m  5  ; 

i-MA  ;  3  =  1;  5  for  cases  which  correspond  to  the  curves  in  Fig. 

11.  At  small  c(  (large  fj  )  function  ^*(1+0*)  depends  weakly  on  c(k 

and  it  is  possible  to  assume  .  At  c(-*oo  function  ¥x(f+fl*) 

approaches  the  limiting  value 

[fft/JP/p0]^  1 +  T+fd«k/4 


§5.  THE  INFLUENCE  OF  HETEROGENEITIES  OF 
FLOW  ON  THE  CHARACTERISTICS  OF  AN  MHD 
GENERATOR  WITH  A  SMALL-SCALE  CHANNEL 


As  an  example  let  us  examine  the  electrical  characteristics  of  an 
MHD  generator  channel  on  the  typical  combustion  products  of  hydro- 

2 

carbon  fuel  [2,  3,  22] .  Let  the  area  of  flow  cross-section  be  'vlO  cm  , 
and  deposits  on  the  walls  of  the  channel  are  absent.  The  relative 
sizes  of  the  channel  and  the  boundary  layers  are  the  following: 

a.  » /;  fi  =3.5;  &3  =  tiu  =  O.W;  ~  =  0.37;  &  *2  /O'* 

Let  us  select  for  evaluation  the  following  parameters  of  gases: 


V:  l  i2  ;  M=2.2;  71 = 2500 °K  ;  P=  2  Catm  (abs.)]j  fe0.5; 

pr  =  0.88 ;  ju  =  5.5- /o's  kg/ms  n  =nr  1  V?  . 

The  induction  of  magnetic  field  B  *  3  ml.  The  temperature  of  the 
walls : 


a)  "cold"  walls  -  fy  *  1000°K;  s  1500°K; 

b)  "hot"  walls  -  T3  ■  2G00°K;  Tu  =  2300°K. 

Figure  13  shows  the  current-voltage  characteristics  for  a  unit  of 
Length  of  channel  at  in  relative  coordinates  y,  y/a^Btu 

at  Z3 *  ta  and  13*X*  (for  - Th/o* °K;  =  f;/o%#  > 

and  in  ^ig.  14  -  the  dependence  of  dimensionless  output  power 
N-N/6  U Z8*Cttf>.  on  current  in  the  load,  calculated  from  the 
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Pig.  13.  Volt-ampere  char< 
acteristics  of  a  model  MHD 

channel  (  fi-7*  *  1000°K, 

Tu  -  1500°K;  5-79  *  2000°K 

Tu  *2300°K). 


Pig.  14.  The  dependence 
of  output  power 


/V~  ttZF&T 

in  the  load 


on  current 

7  J 
*  u-Btt 


us 

1 1 


one-dimensional  method  with  the  utilization  of  a  two— layered  circuit 
of  boundary  layer.  Figure  15  depicts  the  dependence  of  the  short- 
circuiting  current  a  on  the  temperature  of  insulating 

well  at  =  Z0  and  7]  =  1000°K,  1500°K,  2000°K. 


Fig.  15.  The  temperature  effect  of 
insulating  wall  on  the  short-circuiting 
current  „ 

k-aSSbr  (*»"*•)■ 


The  calculations  conducted  showed  that  with  "cold"  channel  iyalls- 
the  output  power  of  the  generator  can  be  2-3  times  less  than  its  “ 
value  not  allowing  for  the  boundary  layer.  The  method  of  calculation1 
of  near-electrode  resistance  exerts  a  strong  influence  on  the  calculated 
value  of  near-electrode  resistance  and  therefore  output  power  at 
T$~  1000°K  (Fig.  1*0.  Since  in  the  given  calculation  the  component 
of  near-electrode  resistance  f}£  was  not  considered,  then  the  results 
obtained  bear  a  semiquantitative  nature.  The  open-circuit  voltage  due 
to  leakages  in  the  turbulent  boundary  layers  decreases  only  by  'v-lO?, 
and  the  short-circuiting  current  falls  by  VL5JS  with  the  reduction  in 
the  temperature  of  the  insulating  walls  from  2000°K  to  1000°K. 
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A  local  circuit  of  substitution  for  an  MHD  generator  channel  has 
been  proposed  which  makes  it  possible  in  a  gas-dynamic  approximation 
to  calculate  the  influence  of  the  heterogeneities  of  flow  on  the 
electrical  characteristics  of  the  generator. 


An  approximate  computation  of  the  basic  elements  of  the  substitu¬ 
tion  circuit  has  been  conducted.  It  has  been  shown  that  in  the  nominal 
regimens  of  operation  of  an  MHD  generator  it  is  possible  to  practically 
disregard  the  leakages  of  current  in  turbulent  boundary  layers  on  the 
insulating  walls.  In  this  case  the  circuit  of  substitution  (Pig.  3) 
is  considerably  simplified  and  the  basic  influence  on  the  character¬ 
istics  of  the  generator  will  be  exerted  only  by  near-electrode  effects, 
especially  the  resistance  of  thermal  boundary  layers,  and  also 
leakage  on  the  material  of  the  insulating  walls  and  the  layer  of 
electro-conductive  deposition  (if  they  are  significant). 


At  low  electrode  temperatures  (  T9  ^  I000°K)  the  selection  of  the 
characteristic  integration  limit  of  t)  exerts  a  powerful  influence  on 
the  magnitude  of  the  components  of  resistance  of  the  boundary  layer 
(f,t)  during  calculation  in  the  two-layered  circuit  of  boundary 
layer,  and  condition  Zj  •  0  can  lead  to  an  increase  of  Rff  by  an 
order.  However,  arbitrariness  in  the  selection  of  Z?  at 
(for  instance,  the  condition  Xt~Zt  at  the  assigned  parameter 


)  in  the  case  %  &  IOOO°R  weakly  influences  the  result 


of  calculation  of  )>  especially  with  the  exponential 


3H 


"-Ji3?  ’*’  V?Ti>?7~ 


approximation  of  the  velocity  (temperature)  profile  up  to  the 
electrode,  and  at  2000CK  it  is  possiole  to  assume  tt  •  0.  This 

apparently  can  explain  the  satisfactory  agreement  of  the  results  of 
the  calculation  of  the  resistance  of  the  thermal  boundary  layer  at 
with  the  experimental  values  obtained  in  work  [9]. 

The  resistance  of  the  thermal  boundary  layer  on  the  electrode  at 
an  assigned  thickness  is  determined  mainly  by  the  temperature  of  the 
inner  surface  of  the  electrode  wall,  especially  at  Z,*Zm  • 

The  structure  of  the  boundary  layer  and  exponent  tl  in  the  case 
of  exponential  approximation  of  the  profiles  of  velocity  and  tempera¬ 
ture  exert  a  noticeable  influence  on  the  real  resistance  of  the 
boundary  layer. 

At  the  large  current  densities  the  profiles  of  velocity  and 
temperature  can  differ  noticeably  from  those  accepted  in  calculation 
and,  furthermore,  a  breakaway  of  electron  temperature  and  the 
emergence  of  a  micro-arc  mode  of  the  flow  of  current  are  possible. 

Thus  the  correct  calculation  of  near-electrode  resistance  requires 
the  careful  analysis  of  processes  which  accompany  the  flow  of  current 
in  the  interelectrode  interval  of  the  channel.  Two-dimensional 
effects  and  near-electrode  resistance  in  the  channel  of  a  Faraday 
generator  with  continuous  electrodes  at  can  cause  a  noticeable 

longitudinal  electric  field,  which  leads  to  an  increase  in  transverse 
conductivity . in  the  interelectrode  interval  of  the  channel.  Evaluations 
show  that  the  Increase  in  output  power  as  compared  with  one-dimensional 
calculation  can  be  significant  (M..5  times  at  fist  1-2,  5). 

The  concrete  calculation  of  a  small-scale  generator  showed  that 
at  an  electrode  temperature  of  ,vl000oK  output  power  can  be  2-3  times 
less  than  its  value  not  allowing  for  the  resistance  of  the  boundary 
layer. 

The  boundary  layer  effect  can  be  significant  also  in  heavy 
installations  at  large  values  of  parameters  Lt/fi  •  If  far  from  the 
separation  point  for  a  hydraulically  smooth  wall  then  at 
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an  identical  number  of  gages  LK/A.  -const  and  an  identical  temperature 

rate  ~  ^  **  ~2 ("7^W)'  .  Since  usually  /77  «s  0.8,  then  the  correlation 

between  the  thickness  of  the  boundary  layer  and  the  transverse 
dimension  of  the  channel  depends  weakly  on  its  length.  Therefore,  for 
the  reduction  of  relative  resistance  in  the  boundary  layer  (or  the 
loss  of  voltage  in  the  channel)  other  conditions  being  equal  it  is 
necessary  to  select  a  relatively  small  number  of  gages. 


The  authors  thank  A.  V.  Gubarev  and  P.  R.  Ulinich  for  useful 
discussions. 
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